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ABSTRACT 

Bright star-forming galaxies radiate well below their Eddington Limit. The value of the flux-mean 
opacity that mediates the radiation force onto matter is orders of magnitude smaller than the UV 
or optical dust opacity. On empirical grounds, it is shown that high-rcdshift ULIRGs radiate at 
two orders of magnitude below their Eddington Limit, while the local starbursters M82 and Arp 220 
radiate at a few percent of their Eddington Limit. A simple model for the radiative transfer of UV 
and optical light in dust-rich environments is considered. Radiation pressure on dust does not greatly 
affect the large-scale gas dynamics of star-forming galaxies. 

Subject headings: galaxies: general galaxies: starburst — dust, extinction — ISM: jets and outflows 
— radiative transfer 



1. INTRODUCTION 

Star formation is inefficient. The timescale at which 
gas is converted into stars is long in comparison to the 
gravitational collapse time. In order to explain this, some 
form of self-regulation or "feedback" is often invoked. 

There are many diflferent mechanisms: core- collapse 
supernovae, thermonuclear supernovae, radiation pres- 
sure on dust, cosmic ray pressure, stellar winds, heating 
due to ionizing radiation, jets from proto-stellar disks, 
turbulent stresses, magnetic stresses, ram-pressure strip- 
ping, pulsars winds, cosmic ray heating, dark matter an- 
nihilation as well as the various forms of energy release 
from super-massive black hole growth (see the Appendix 
for a brief literature review) . 

The UV and optical radiation force tliat directly re- 
sults from star formation seems to be a promising can- 
didate and has recently gained a significant amount of 
attention. Second only to prompt neutrinos from core- 
collapse supernovae, starlight is the dominant source of 
energy release in the Universe. The nuclear binding en- 
ergy that efficiently fuels starlight is orders of magni- 
tude larger than the gravitational binding energy of even 
the most massive galaxies. In addition, UV and opti- 
cal starlight is extremely well coupled to the interstellar 
medium via scattering and absorption onto dust grains. 
An assessment of this mechanism is the subject of what 
follows. 

2. THE EDDINGTON LIMIT 

The radiation force density is given by (cf. Blaes & 
Socrates 2003; Rybicki & Lightman 1979) 



/ 



du Ku ¥i, 



(1) 



where is the the sum of scattering and absorption 
opacities at photon frequency v, Fj^ is the radiative flux 

and p is the fluid mass density. By equating this to the 
gravitational force and integrating over the galactic sur- 
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face we arrive at the Eddington Limit 
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(2) 



where Menc is the enclosed dynamical mass and is the 
flux-mean opacity, which can be defined as 



J dv 
J duLi, ' 



(3) 



if the Spectral Energy Distribution (SED) and the opac- 
ity law are uniform on the surface of the galaxy, so that 



J dA ■ j du Kv^v = j j 



dA- K^F^ = L . (4) 



Note that the Eddington Limit is a statement of hy- 
drostatic balance that is independent of geometry and 
whether or not the flow is optically thick or thin (Socrates 
2012; Abramowicz et al. 1980). 

The most luminous starburst galaxies in the Universe 
have luminosities of order Lmax ^ 10^^ ergs/s. The 
starlight couples to gas by scattering and absorption onto 
dust grains. The UV opacity on dust is k^^ ~ lO^Kes, 
where Kgs is the electron scattering opacity. The most 
massive galaxies in the Universe possess an enclosed mass 
of order Mr.,-,c 



10^^M(^. With this combination of 



L 



max J /^uv ^^'^ Mono, it sccms as though L„ 



K^^. This serves as the central ar- 



when choosing 

gument behind the Murray et al. (2005) photon feedback 
model. 

Such an argument is appealing. It provides a simple 
and compelling theoretical framework for understanding 
the coupling of the various forms of energy release that 
accompany star formation with the gravitational content 
of the galaxy in question. Ultimately, many of the prin- 
ciples outlined in Murray et al. (2005) may contribute 
to our understanding of the Faber-Jackson (1976) rela- 
tion, the fundamental plane of elliptical galaxies, and the 
maximum luminosity of star-forming galaxies. 

However, k^^ Kf- f^-ct, for starbursting galaxies 
Kp <C Kjjy. In Figure 1, the frequency, or wavelength, 
dependent dust opacity (solid red line) is super-imposed 
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upon the SED of the compact ~ presumably - starburst- 
ing galaxies Arp 220 and M82, and upon an average spec- 
trum of high-z ULIRGs (black solid line) . The flux- mean 
opacity ^ 5 — 20cm^/g for all of the sources (red 
dashed line). Bright star-forming galaxies radiate well 
below their photon Eddington Limit (blue dotted line). 

The primary cause of this shortfall in the radiation 
force is that most of the light belongs to wavelengths 
where the opacity is small. Or conversely, there is rela- 
tively little light at frequencies where the opacity is large. 

The situation is similar to the radiative transfer in a 
stellar envelope where the opacity is large for values of 
photon energy near atomic transitions. Yet, photons in 
stellar envelopes escape primarily by avoiding such re- 
gions of high opacity, which leads to deep absorption 
features. Consequently, the flux-mean opacity in stellar 
envelopes remains close to the low continuum value. 

3. UV/OPTICAL SCREENING 

Figure 1 indicates that the ratio of the UV to FIR lu- 
minosity in bright star-forming galaxies is low - of order 
a few percent. As a result, the flux- mean opacity is 
correspondingly diminished and the galaxy in question 
radiates well below its Eddington value. 

An increase in the gas surface density leads to a corre- 
sponding increase in the rate of star formation. There- 
fore, dense gas serves both as a source of UV/optical 
starlight radiation pressure via star formation, while 
simultaneously screening itself via absorption and re- 
processing of the UV light ultimately into the FIR. A 
depiction of the UV/optical screening in a dense star- 
forming galaxy is shown in Figure 2. 

Consider a galaxy of characteristic length scale i?, with 
uniform density, temperature (p, T) and composition. 
Under these assumptions, the starlight screening length 
/ is a constant. For a constant value of the star forma- 
tion rate per unit volume, the ratio of UV light to FIR 
radiation is 

^ ~ % ■ (5) 

The expression above may be a good approximation for 
cool dust-rich galaxies that are opaque to starlight, where 
the FIR luminosity Lpjj^ is responsible for most of the 
bolometric energy release. 

The simple picture outlined above may explain the 
small ratio of UV to FIR photon power in starburst- 
ing galaxies. An attempt to quantify the properties of 
starlight screening by considering the equation of radia- 
tive transfer is given below. 

3.1. Radiative Transfer 

The spatial and directional dependence of the starlight 
is quantified by the specific intensity I^. Dust grains 
are the primary source of scattering and absorption of 
starlight, where absorption acts to redistribute the input 
stellar photon energy longwards to IR and FIR frequen- 
cies. At these wavelengths, the dust opacity is relatively 
small and the radiant energy can more easily escape. 

In steady state, the equation of radiative transfer for 
the specific intensity I^ix^s) along a direction s reads 
(Rybicki and Lightman 1979) 

V • [Iu{s)s\ = -{a^ + a^)I^{S) + 
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Fig. 1. — SED (black solid line) and total luminosity (black 
dashed line), wavelength-dependent opacity kx (red solid line, in- 
cluding both scattering and absorption) and flux-mean opacity (red 
dashed line) and the Eddington Limit (blue dotted line) for lumi- 
nous star-forming galaxies. For M82 and Arp 220, the SED is 
taken from Silva et al. (1998). The enclosed dynamical mass is 
estimated as Mcnc ~ 2 X IO^Mq for M82 (Greco et al. 2012) and 
A/cnc ~ 4 X 10"' M0 for Arp 220 (Silva et al. 1998). The average 
SED and median enclosed mass for ULIRGs are taken from the 
sample of Michalowski et al. (2010). The opacity law is assumed 
to follow that of the Milky Way as calculated by Weingartner & 
Draine (2001) for Ry = 3.1 and a gas-to-dust mass ratio fg^ = 100. 
The opacity model by Semenov et al. (2003) yields similar results. 



+<T,J Us')u{s-s')dn' + J,{S) (6) 

where the spatial dependence on x has been suppressed. 
Here, a^{x) is the absorption coefficient and a,^{x) the 
scattering coefficient. They are, respectively, equal to 
the inverse of the mean free paths for absorption and 
scattering of photons of frequency v. The emission coef- 
ficient, or source, is given by ju{x, s) and the normalized 
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screening 
length 



Fig. 2. — A simple model of a star-forming galaxy. The starlight 
is assumed to be injected uniformly within the radius R. Due to 
the large optical depth at UV/optical wavelengths, the cold gas 
in the bulk of the galaxy screens the starlight except for an outer 
region whose thickness is of order the screening length I. The bolo- 
metric output generated from within the galaxy ultimately escapes 
in the form of quasi-thermal FIR radiation. Similar arguments for 
explaining the ratio of UV/FIR power can be made in cylindrical 
geometry as well, with the same result. 



differential scattering cross section, /^(s • s), satisfies 



/ f^{s- s')dn' = 1 



(7) 



By defining the mean intensity Ji,{x) = J Iv{x,s) d^, 
the radiative flux Y^{x) — J Iiy{x,s) sdi} and the mean 

emissivity q,y{x) = J j^{x, s) dfl,^ the radiative transfer 
equation averaged over the propagation angle reduces to 



V-F, 



(8) 



In the diffusion approximation, the specific intensity I^, 
can be expressed as an isotropic part (i.e., the mean in- 
tensity Jv) plus a small directional flux F,^. That is, 



I,{s) = J, + -^F, 
47r 47r 



(9) 



where the factor of three on the right hand side re- 
sults from enforcing the Eddington closure approxima- 
tion. Substitution of this expression into eq. (6), multi- 
plying by s and integrating over all solid angles, leads to 
a relation between the mean intensity and the flux of the 
following form 



F^ = -D^VJ^ 



where 



(10) 



(11) 



is the photon diffusion coefficient, and g is the average co- 
sine of the scattering angle. The two relations in eqs. (8) 
and (10) between the mean intensity and the flux can 

^ In the following, the source, quantified by the emission coeffi- 
cient jv , is taken to emit isotropically. 



be combined to give a diffusion equation for the mean 
intensity alone 



V • (£)^V J^) - a„J„ 



(12) 



In a uniform medium, the diffusion coefficient Di, is 
independent of position, and under these conditions 



-qi' 



(13) 



which is in the form of a screened Poisson's equation 
with screening length l^, = sjl\ja^. The Green's func- 
tion Gy[x,x'^ for the screened Poisson's equation satis- 
fies V^Gi/ — Gu/ll = —6{x — x') and is given by 



G„{x,x') 



exp(-|a; - x'\/l,y) 



Awlx ■ 



X' 



(14) 



which is also referred to as the Yukawa potential. Un- 
der the assumptions that the radiation field is quasi- 
isotropic and photon diffusion proceeds through a uni- 
form medium, the analogy of the radiative transfer prob- 
lem with electrostatics is clear. The average intensity 
Jj/, which is equivalent to the photon energy density, can 
be thought of as an electrostatic potential that is pro- 
duced by some source of radiation q„, which itself can be 
thought of as an electrostatic charge density. Further- 
more, it follows that the radiative fiux Fj^, responsible 
for the radiation force, is analogous to an electric field. 
Absorption attenuates the mean intensity in starlight, 
which is therefore evanescent in the bulk of the galaxy, 
and its effect may be thought of as a screening. 

The solution to the radiative transfer of starlight un- 
der these circumstances further motivates the arguments 
behind Figure 2 and eq. (5). 

3.2. Boundary Value Problem in Spherical Symmetry 

In a uniform sphere, the source function = qo is 
a constant. The inhomogeneous solution for J^, to the 
screened diffusion equation must be regular at the center 
of the galaxy and satisfy an appropriate boundary con- 
dition at the surface r = R. For the latter, we employ 
the two-stream approximation (Rybicki and Lightman 
1979), which relates the flux F^ to the mean intensity 
in such a way as to maintain the Eddington approx- 
imation for the starlight. With this, the mean intensity 
inside the spherical galaxy is given by 



go il 



1 - 



R smh{r/li,) 

7 c 



(15) 



where the constant C = C{R) is determined by the 
boundary condition at r = i? and is written as 



C=^/3D^cosh{ ^ ]+{ 1-V3 



R 



sinh ( ^ 



.(16) 



For li, <C R, the mean intensity inside the sphere is al- 
most constant and equal to ~ q^tl/Di, ^ q^^/av, as 
expected for a homogeneous source with emissivity go 
that is absorbed locally with absorption coefficient a^,. 

The radial component of the radiative flux = r • F^ 
is given by 



qoRlv 



r cosh(r/Zy) — li, sinh(r//i^) 



C 



, (17) 
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which vanishes in the hmit r — > 0, as expected. Near 
r = 0, iv rises hnearly with r up to r « 1^, at whicli 
point Fi, increases, approximately, exponentially with 
scale length li, up until r = R, where it reaches its 
maximum value. Since the UV radiation force is pro- 
portional to Fj,, it follows that, in the limit l^^/R « 
L^jy/Lpjp <C 1, throughout the bulk of the flow the UV ra- 
diation force, effectively, vanishes. Conscqiiently, the so- 
called "single-scattering" approximation for UV photons 
(Murray, Quataert & Thompson 2005; 2010; cf. Hopkins 
et al. 2012 a,b,c) significantly - or exponentially - over- 
estimates the UV radiation force throughout the bulk of 
a given star forming flow. 

In the absence of absorption (set a,, = in eq. (13)), 
while still enforcing the two-stream approximation at the 

surface, the un-screened intensity ji"' is given by 



go 



(i?2 



qoR 
^/3 



(18) 



where D^^ — [3(1 — g)(T^]~^ ^ 2_D^ for isotropic diffusion 
(i.e., g = 0) and assuming equal mean free paths for scat- 
tering and absorption (i.e., cr^ = a^). The corresponding 

un-screened radiative flux is then F^^^ = qor/3. 

For R ^ and R ^ Di, (which also implies R ^ 
l^), the ratio between screened and un-screened mean 
intensities at the center (i.e., at r = 0) is 



7(0) 



ll 



r=0 



12/2 



R^ i?2 



(19) 



It follows that, for galaxies such as Arp 220 and ULIRGs 
whose starlight is highly screened, the UV/ optical radia- 
tion pressure is strongly suppressed at depth. 

The ratio of outward radial fluxes at the outer bound- 
ary of the galaxy (i.e., &t r = R) reads 



(0) 



1 /j/ 



(20) 



which corroborates the intuition that led to Figure 2. 



Galaxies, such as the Milky Way, that are not as rich in 
cool dense molecular gas do not suffer from the same level 
of self-shielding that renders radiation pressure ineffec- 
tive in starbursting sources. However, the lack of molec- 
ular fuel also implies a relatively low level of star forma- 
tion. Although the flux-mean opacity may be higher, due 
to a harder spectrum, the total luminosity is low. So, it 
is unlikely that un-screened galaxies radiate above their 
Eddington Limit. 

The radiative transfer of starlight in luminous star- 
forming galaxies can be modeled in the diffusion ap- 
proximation, including the effect of absorption by dust 
grains. Dust-screening greatly diminishes the value of the 
starlight intensity, or pressure, at depth. Only an outer 
shell with thickness equal to the screening length con- 
tributes to the outgoing flux of UV/optical starlight. An 
increase in the density of molecular gas increases the rate 
of starlight production and the Eddington Limit may be 
approached. Yet. increasing the fuel supply necessarily 
leads to an increase in dust-screening, with correspond- 
ing downward departures from the Eddington Limit, due 
to the corresponding decrease in the flux-mean opacity. 

Since the sources under consideration radiate well be- 
low their Eddington Limit, the large-scale radiation force 
due to FIR photons is also small in comparison to the 
galaxy-scale gravitational force at the surface. It is there- 
fore unlikely that the FIR radiation pressure will be large 
enough to be dynamically important at depth. In the dif- 
fusion approximation, radiation pressure is proportional 
to optical depth, which is proportional to column density. 
Since gas pressure is also proportional to column density, 
the ratio between gas and FIR radiation pressure should 
remain approximately a constant with depth. 

There may be some class of galaxies in the Universe or 
objects within luminous star-forming galaxies themselves 
(e.g., giant molecular clouds), where starlight pressure on 
interstellar dust leads to a radiation force that exceeds 
the local source of gravity (Murray, Quataert & Thomp- 
son 2010). It may be useful to empirically determine 
the Eddington Limit for these objects with the method 
described in §2. 



4. DISCUSSION AND SUMMARY 

In §2 an empirical method - the calculation of the Ed- 
dington Limit - for determining the relative strength of 
the radiation force on the surface of galaxies is described. 
It is demonstrated that light resulting from star forma- 
tion is too dim to affect the large-scale gas dynamics of 
actively star-forming galaxies. 

The dust opacity is large at small wavelengths, where 
the emergent spectrum is relatively dark, while the dust 
opacity is small at the long photon wavelengths, where 
the emergent spectrum is bright. By combining these 
two properties it is apparent that the radiation force is, 
in general, unexpectedly weak in gas-rich star-forming 
sources. 
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APPENDIX 
REVIEW OF RADIATION FEEDBACK 

In recent years, the idea that radiation pressure due to starUght interacting with dust grains may be comparable to 
gravity was considered by Scoville (2003). Martin (2004) and Murray et al. (2005; MQT) hypothesized that radiation 
pressure on dust grains may be responsible for driving large-scale galactic winds in analogy to line-driven stellar winds 
of massive stars. MQT utilize the large UV/optical opacity on dust to justify that, despite the low light to mass ratio 
of the galaxies, it is possible for them to radiate with super-Eddington luminosities. Thompson et al. (2005) focus 
on the effects of coupling the IR and FIR radiation to dust and thus, to the gaseous component, in an attempt to 
construct a starlight-powered galactic disk solution. More recently, Murray et al. (2010) changed focus from large 
scales to the scale of giant molecular clouds (GMCs), where they posit that fully-populated star clusters born in GMCs 
unbind the site of their birth once the star to gas ratio becomes large. In a series of recent papers (Hopkins et al. 2011; 
Hopkins et al. 2012a, b,c; Faucher-Giguere et al. 2013), both the GMC-scale starlight-driven hypothesis as well as an 
approximate formulation of the large-scale radiation force are utilized in a number of massive numerical simulations 
meant to describe a variety of star-forming galaxies. They find that the radiation force is capable of significantly 
altering the gas dynamics of star-forming galaxies. 

Recently, Andrews & Thompson (2011) and Zhang & Thompson (2012) note that, at face value, star-forming galaxies 
are sub-Eddington with respect to the Rosseland mean opacity for dust. By taking the Rosseland mean, they assume 
a functional form for the SED of the galaxy. As a result of this, they are led to believe that there is some level 
of theoretical ambiguity when assessing the strength of the radiative forcing, particularly when discerning between 
optically thick and optically thin regimes. One of the most attractive features of the Eddington limit is that it is an 
empirical tool. In §2 the observed SED of galaxies is integrated against a typical dust opacity law in order to obtain 
the strength of the radiative forcing on the surface of a given galaxy, with relatively little ambiguity. Furthermore, the 
Eddington Limit is independent of whether or not the flow is optically thick or thin. 

In describing the feedback of cosmic ray protons during star formation, Socrates et al. (2008; see also Sironi & 
Socrates 2010 for the case of quasars) attempt to draw a distinction between cosmic ray feedback and the feedback 
from starlight. In their appendix, they note that UV photons decouple and are reprocessed into the FIR close to the 
massive star of their origin. Once starlight is converted into the FIR band, they argue that it is unlikely to affect 
the dynamics of the interstellar medium. A critical re-examination of photon feedback is the subject of many recent 
numerical studies (Novak et al. 2012; Krumholz & Thomson 2012; Wise et al. 2012; Agertz et al. 2012; Jiang et 
al. 2013). They all concur, to varying degrees and for various reasons, with Socrates et al. (2008) in that that 
starlight pressure may not be efficient in its coupling to the interstellar medium. However, to solve the full radiation 
hydrodynamical problem of photon feedback is incredibly challenging from a numerical perspective due to large ratio 
of the photon diffusion to gas sound speed (cf. Fernandez & Socrates 2012). In §2 an attempt is made to circumvent 
this debate by empirically determining the strength of the radiative forcing on the surface of star-forming galaxies. 



